On Feynman's calculation of the Frohlich polaron mass 
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Feynman's formula for the effective mass of the Frohlich polaron is rederived from the formalism 
of projected partition functions. The mass is calculated as inverse of the diffusion coefficient of the 
polaron trajectory in imaginary time. It is shown that correlation between the electron and phonon 
boundary conditions in imaginary time is necessary for consistent derivation of the Feynman result. 



PACS numbers: 71.38.Fp 



INTRODUCTION 



Fifty years ago Feynman published his seminal paperi 
on the Frohlich polaron, in which he combined path inte- 
gration and an action variational principle to obtain the 
polaron energy and effective mass. Subsequently, Feyn- 
man's method was generalized by many authorsS^*^, 
for comprehensive review see for example Refs.[H||. 
When the exact polaron energy and mass were calculated 
by Fourier Quantum Monte Carlc^IS and Diagrammatic 
Quantum Monte CarloAi*i£ methods, Feynman's polaron 
energy was found to be remarkably accurate deviating 
from the exact value by less than 0.5% for all couplings. 
In contrast, Feynman's mass formulaic overestimated 
the exact mass by as much as 50% leaving room for 
possible improvements. Generalizations of the original 
calculationSiii^ii^ produced either no or very small, less 
than 1%, numerical correction to the mass. (In some in- 
stances, the correction was of the wrong sign&ii) Thus 
the generalizations did not improve the agreement be- 
tween the analytical and numerical masses. The much 
larger error in the mass was not regarded as something 
unexpected. It was known since Feynman's paper that 
the mass did not satisfy a variational principle, and there- 
fore the accuracy of the approximate treatment was ex- 
pected to be somewhat uncontrolled. 

In Feynman's methocU*A£iil, phonon variables are in- 
tegrated out analytically resulting in a self-interacting 
retarded one-electron action. Phonon integration is per- 
formed under periodic boundary conditions in imaginary 
time. Periodic boundary conditions are also assumed for 
the electron trajectory. Both conditions are perfectly 
appropriate for the energy calculation, since the ground 
state energy is obtained from a full thermodynamic par- 
tition function. However, later, in the mass calculation 
the electron trajectory is broken and its r = and r = 
ends are displaced relative to each other. The effective 
mass is inferred from variation of the polaron energy with 
real space separation of the two ends of the trajectory. 
In Feynman's approach, the retarded polaron action is 
not corrected for this displacement. In other words, the 
action resulted from phonon integration for a periodic 
electron trajectory is carried over to an open electron 
trajectory. 



was realized that the boundary conditions for the phonon 
and electron variables must be synchronized for consis- 
tent calculation of the polaron mass. That is, they are 
either both periodic or both open. And when they are 
both open they are correlated. The correlation follows 
because the electron and phonons are coupled in one sys- 
tem and share a common integral of motion, the polaron 
momentum. Accurate account of the correlation results 
in a polaron action different from the periodic phonon 
action derived by Feynman. Thus the question arises of 
whether this difference is significant enough to yield a 
correction to the Feynman mass formula. Such a cor- 
rection might account for the deviation from the Monte 
Carlo mass. 

This question is investigated in the present paper. It 
will be shown that the correlation between the phonon 
and electron boundary conditions is essential for consis- 
tency of the mass calculation. Without it, a second dou- 
ble integral arising from phonon integration cannot be ne- 
glected in the mass calculation while it can be neglected 
in the energy case. Moreover, it spoils the entire mass 
calculation. The correlation between the two boundary 
conditions restores internal consistency and leads to the 
correct final formula. It will also be shown that these ob- 
servations produce no numerical correction to the Feyn- 
man result. 



II. EFFECTIVE MASS FROM SHIFTED 
BOUNDARY CONDITIONS 

The polaron mass can be calculated as the inverse dif- 
fusion coefficient of an open-ended polaron path in imag- 
inary time. This section contains a derivation of this 
relation. The derivation is general and valid for any non- 
relativistic composite particle. The momentum KK. of a 
translation invariant system is a constant of motion. In 
the polaron case, this is the sum of an electron momen- 
tum and momenta of all excited optical phonons. One 
defines the projected partition function as a Gibbs sum 
restricted to states with the same K: 
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Here |n) are eigenstates of the Hamiltonian H, KK n is 



it the momentum of state |n), and f3 — (fc^T) is inverse 
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temperature. The system is assumed to occupy a finite 
volume V, and the wave vector is quantized to a discrete 
set of values K„. To transform Zjc, introduce real-space 
configurations |Q) which are direct products of all the de- 
grees of freedom in the system. For the Frohlich polaron, 
\Q) = l r ) Ilr' l-P( r '))i where r is the electron position 
and P(r') is the polarization at point r'. The states Q 
form a complete orthogonal basis, I — J dQ\Q)(Q\, and 
(Q1IQ2) = &{Qi — Q2), where the unity operator and 
the delta function are also direct products between dif- 
ferent degrees of freedom. Inserting two unity operators 
in Eq. Q the partition function is rewritten as follows 



dQ 1 dQ 2 {Q 2 \e-' 3H \Qi)-WK 



W K = 5^<Qi|n)<n|Q 2 ) <5kk„ = (Qi|K)(K|Q 2 



(2) 
(3) 



The meaning of the last expression is that both con- 
figurations Qi and Q 2 have to be projected on the 
same wave vector K. To perform projection consider 
a parallel shift of configuration Q by a three dimen- 
sional vector R. The resulting state will be denoted 
as \Q + R). (Note that summation is only symbolic 
here.) An arbitrary configuration Q generates a fam- 
ily of states \K Q ) = V- 1 J dR e" iKR |Q + R). Inversely, 
\Q + R) = £ K e iKR |K Q ). In Wk only the respective 
components IKqJ and |Kq 2 ) of \Q\) and \Q%) survive 
projection on K. As a result, one obtains 
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K 



(K 0l IK, 



= ^ J dR 1 dR 2 {Q 1 +R 1 |Q 2 + R 2 )e 4K ( Rl - R2 ) 

= i J d(Ar)e lKAr (Q 1 + Ar|Q 2 ) 

= i J d(ArK KAr • 5 [{Q x + Ar) - Q 2 ] , (4) 

where Ar = Ri — R 2 . The resulting delta function en- 
sures that a many body configuration Q 2 is identical to a 
configuration Q\ shifted by Ar. Substitution in Eq. 
and integration over Q 2 yield 



Z ^ = \\ d(Av)e lKAr J dQ(Q + Ai 



-0H 



\Q)- (5) 



The matrix element under the dQ integral is the density 
matrix operator taken between an arbitrary real-space 
configuration Q and the same configuration shifted by 
Ar. Since Q is a many body state all particles shift in 
parallel. In the polaron system the electron coordinate 
r and polarization profile P(r') shift together. Equation 
© suggests defining the shift partition function 



dQ{Q + Ar\e- 0H \Q) 



(6) 



which is characterized by the shift vector Ar. The zero- 
shift partition function, ^Ar=0j coincides with the usual 



thermodynamic partition function. Equation J5J states 
that the projected partition function and the shift parti- 
tion function satisfy a Fourier-type relation 
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(7) 



This relation is valid at any temperature T. In the low 
temperature limit, Zk is dominated by the lowest energy 
eigenvalue with given wave vector, which allows deriva- 
tion of a useful formula for the effective mass. At small 
K, the energy is approximated by £k = Eq + , and 
the projected partition function Zyl — > e~° Eli as T — > 0. 
Expanding Eq. Q to the second order in K one obtains 
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Z A r-(8) 



On the right, the linear term in K vanishes after integra- 
tion by inversion symmetry: Z-Ar = ZAr- The rest is 
transformed as follows 



/3ft 2 K 2 _ MAr)(KAr) 2 Z Ar = 2 , 
m* fd(Ar)Z Al ~ U 1 ■ 



shift ■ 



(9) 



The definition of Zat, Eq. iJBJ, implies that the ratio of 
two integrals in the last expression is the mean value of 
(KAr) 2 evaluated with shifted boundary conditions in 
imaginary time. The latter means the initial (at imagi- 
nary time = 0) and final (at imaginary time = (3) con- 
figurations are the same, cf. Eq. JfjJ, but they can be 
shifted relative to each other by a three dimensional vec- 
tor Ar. This shift vector is arbitrary. Averaging under 
shifted boundary conditions is understood hereafter as 
averaging over Ar, with the weight given by Zav Upon 
expansion of the square in Eq. @ , the mixed terms aver- 
age to zero by symmetry, ((Arj)(Arj)) S hift = 0, while the 
diagonal terms are equal, ((Ar l ) 2 ) shi f t = i((Ar) 2 ) shift . 



That results in 



1 _ 1 /d(Ar)(Ar) 2 Z Ar _ ((Ar) 2 ) shift 



3(3h 2 J d{Ar)Z Ar 



3(3h 2 



(10) 



This equation allows an elegant interpretation of the ef- 
fective mass in terms of imaginary time diffusion. Since 
the shift vector is not fixed, the system evolution from the 
initial to the final configuration can be regarded as dif- 
fusion during time t = hfj. In normal three-dimensional 
diffusion, the mean squared displacement is proportional 
to the time interval, ((Ar) 2 ) = 6Dt, where D is the dif- 
fusion coefficient. Thus Eq. I|10|l is rewritten as 
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m* 



2D 



(11) 



Note, that Eq. I|10|l can also be regarded as a fluctuation- 
dissipation relation. The effective mass characterizes dy- 
namical response, while the mean squared displacement 
is an equilibrium property. 
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Equation jTDj) is especially useful in understanding 
mass enhancement of composite particles such as the po- 
laron. Interaction with various fields (e.g., phonons) in- 
creases the statistical weight of trajectories with small 
((Ar) 2 ) thereby slowing down the diffusion and increas- 
ing the particle's mass. Thus the mass enhancement is 
conveniently visualized as the increased "stiffness" of the 
trajectories. Equation (jTU|l or its analogues were used 
in Monte Carlo calculations of the effective masses of 
polarons&ifliiSiiiiSi, bipolaronsSi, and defects in super- 
fluid helium 2 ^. 



III. 



POLARON ACTION 



The results of the previous section have important im- 
plications for the Frohlich polaron. As soon as the two 
ends of the electron path are allowed to shift relative 
from each other to obtain the mass, the polarization pro- 
file must shift accordingly. Thus phonon integration has 
to be performed under more general boundary conditions 
in imaginary time than periodic. In general, this should 
modify the polaron action. The modified polaron action 
is calculated in this section. 

A starting point is the polaron action as formulated by 
Frohlich 23 and Schultaife 
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S[r(r);P(r',r)] =- / dr 
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(12) 



(13) 



Here 51 is the optical phonon frequency, Eq and £oo are 
the static and high-frequency electric permittivities of the 



crystal, |e| is the unit charge, and m is the band mass of 
the electron. r(r) is the imaginary-time electron trajec- 
tory. P(r',r) is the imaginary-time polarization trajec- 
tory. A dot above a variable denotes partial derivative 
with respect to imaginary time r. The action is supple- 
mented by shifted boundary conditions 



r(/3) = r(0) + Ar, 
P(r',/3) = P(r'-Ar,0), 



(14) 
(15) 



where Ar is an arbitrary 3-dimensional vector. Except 
for the boundary conditions phonon integration proceeds 
along the lines outlined by Schultz 16 . The polarization 
field is expanded in a Fourier series with real amplitudes 
A and B : 

P(r',r) = M T^^ q (r)cosqr' + S q (r)sinqr']. 

( fe >o) |q| 

Note that the sum over q extends only over half of mo- 
mentum space, which is indicated by l (q x > 0)'. The 
Fourier-transformed polaron action 1|12[) and boundary 
condition (|15fl become 



5[r(r); J 4 q (r), J B q (r)] = - / dr— v 

Jo *n 
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+ IqT V I [Aq(T) qr(T) " Sq(T) C ° S qr(T) T ' 



A q (/3)cosqAr + B q (/3) sinqAr = A q (0) (18) 
-A q (/3)sinqAr + B q (/3) cosqAr = B q (0). (19) 

Since action l|17|) is diagonal in amplitudes A and B, 
path integration can be performed for each component 
independently. Using the standard methodsiii&ii one 
obtains in the low-temperature limit e hni3 ^> 1: 



J 



S[r(r); J 4 q (0),A q ( / 3), J B q (0), J B q (/3)] = - / dr^- + £ S q , 
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mr 



(20) 



Sq = --^M [Al(0) + Al(P) + Bl(0) + Bl((3)] +!L&!WL1JJ o dr'rfr"G(r',r")cosq[r(r')-r(r")] + 

+ (47rH] JT (_ mr [A ^ Q) sinqr(r) _ Bq(Q) cosqr(r)] + e -MJ03-r) ^ q(/3) sinqr(r ) _ Bq(/3 ) cosqr(r)]} , (21) 
q V VJq 1 J 



G(r',r") 



sinhMV ■ sinh hn((3- r"); r' < t" 

Ml sinh hto0 ' \ sinh7ifi(/3 - t') ■ sinhMV'; r' > r" 



(22) 



The action is still a functional of two end polarizations, at r = and r = (3. However, those are related by the 
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conditions l|18|) - (|19|) . Final integration over the end variables leads, after straightforward algebra, to 



SavW)] = dT ^2+ TTTt; / / drdT 



ii 



2Ti 2 2V2\mJJ J ) |r(r')-r(r")| |[r(r') - r(r")] sgn(r' - r") - Ar| 



1 1 \ 2mn 



2TiQ, \eoo e 



(23) 



(24) 



where a is the Frohlich coupling constant. Ar in the denominator of the last term in Eq. is a direct consequence 
of the shifted boundary conditions for the polarization field. It is convenient to perform a linear transformation 
r(r) = r'(r) + ^Ar, which transforms path integration to periodic boundary conditions. The resulting action is 



m (Ar) 2 

2ft 2- 3~ Jo 2h 2 



Fi r / \i >"> y-* 1 J / i 1111 

SAr[r(r)]=--5^-- / dr- J + 



_ W ff dT < dT n J e -Rn|r'-r"| e _ fi n(^-|r'-r"|) 1 

+ 2V2V / ioio T T \|r(r')-r(r") + ^lAr| + |[ r ( T /) _ r ( r ")] sgn ( r ' - T /') - £^=I^Ar| J ' 



r 



The bar over 5* indicates that the action is a functional of 
a path periodic in imaginary time. The full shift partition 
function JfjJ) is given by 

r(r,0) 

Z A r = Z ph - Vve s ^ r ^\ (26) 

where T)v is path integration over electron coordinates. 
.Zph is the partition function of a free polarization field. 
This is a multiplicative constant that cancels out in the 
mass calculation, cf. Eq. l]lUp. 

The polaron mass originates from explicit dependence 
of action (|25|l on the shift vector Ar. The first term 
corresponds to the bare electron mass. Phonon-induced 
mass enhancement comes from the double integral. The 
two integrands have similar functional dependence on 
r(r), but make different contributions to the action. The 
first integrand exponentially decays away from the di- 
agonal r' = t" . Since the odd powers of Ar vanish 
in path integration, the first fraction's contribution is 
O{0) +0((Ar) 2 p- 1 ) + 0((Ar) 4 p- 3 ) + ... The first term 
of this expansion adds to the polaron energy, while the 
second one adds to the mass. 

Consider the second fraction in the double integral 
(I25|l . which is the main focus of the present study. The 
exponential numerator limits integration to finite inter- 
vals around the points (0, 0) and (/3,0). In the denom- 
inator, the combination (/3 — |t' — t"\) = 0(1). There- 
fore the second fraction's contribution to the action is 
0(1) + 0((Ar) 2 /T 2 ) +0((Ar) 4 /?- 4 ) + . . . In the /3 -► oo 
limit, each term in this expansion is small in compari- 
son with the corresponding term from the previous series 
with the same power of Ar. Thus the entire second frac- 
tion can be omitted in favor of the first one. 

Here comes the critical observation. Such a nice term- 
by-term domination of the first fraction over the second 



one takes place only as a result of the shifted boundary 
conditions for the polarization. Indeed, without the lat- 
ter the combination in the denominator of the second 
fraction in Eq. 125|) would have been ^ T ~ T I Ar = 0(Ar) 

instead of it±±^\) Ar = ©((Ar)/?" 1 ). As a result, the 
second expansion would have been 0(1) + 0((Ar) 2 ) + 
0((Ar) 4 ) + . . . The first term is still small compared 
to the corresponding term from the first fraction, which 
leads to the correct polaron energy. However, the (Ar) 2 
term now dominates its counterpart from the first frac- 
tion, which totally confuses calculation of the effective 
mass. Thus the neglect of the shifted boundary condi- 
tions in phonon integration results in a serious internal 
inconsistency in the mass calculation. Apparently, Feyn- 
man avoided this difficulty by omitting the second frac- 
tion in Eq. (|23|l from the outset. Had he retained the 
full form of the phonon propagator, including the second 
part e~ riS1 ^~' T ~ r ^ , he would have faced the problem 
outlined here. 

The analysis presented in this section enables to safely 
neglect the second fraction in Eq. (|25l) . However, it will 
still be included in the forthcoming mass calculation in 
order to illustrate further the above argument. 



IV. POLARON MASS 

Polaron action (|25() is real. Therefore the shifted par- 
tition function satisfies the Jensen-Feynman inequality 

r° Vj . . e S Ar [r(r)] > e <5 Ar [r(T)]-SL[r(T)]} . 

. r 'a-.e^l'M], (27) 

J r 
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where (...)o denotes averaging with the trial action 5*0. 
Feynman's trial model consists of two particles with 
masses m and M, which are elastically coupled with 
a spring constant k. Note that if the first particle's 
mass is different from the electron mass m, the differ- 
ence (S — §o)o diverges. The second particle's mass 
M and k are variational parameters. It is customary 
to replace them with two parameters w — ^\f^ and 

v = w-\/l + ^- which both have units of energy. Calcu- 



lation of the right-hand side of Eq. (|27|l is tedious, but 
proceeds along essentially the same lines as the original 
Feynman calculation. Therefore the intermediate steps 
are not presented here. An important note concerns the 
Ar dependence of the trial action. It is fully represented 
by the term — 2S^"^ir - > reflecting the fact that the 
total mass of the trial model is m + M . As a result of 
the calculation, inequality i|27|) takes the form 



/ m 



„x P { —(3Ep + 0(1) - + 0((Ar) 4 /T 3 ) 



(28) 



E F = - 



3 (v — w) 2 



4 v 



l [P [P e -nn\r'-T"\ _|_ e -m(.p-\T'— r"\) 

[*(r',r")]' ' 



— = (7ift)S~ / / dr'dT 1 ' 



(29) 



a si fPf 13 \ T ' _ T n\2 e -hO.\T'-T"\ , fg_l /_ /m2 -hn(0-\r'-T"\) 1 

m F = m{ l + J?L(ftfi)ii / / dr'dr"^ ^ + lP |T 3 T U I, (30) 



6^ " /? 7 7o 



[$( r ' ;T ")]T 



$(r',r") = 



A/ 



(r' - t") 2 



/5 



1 _ I e ~ 2l,r ' - i e -2"(/5-r') 



2v(/3-t") + e _„( T ' +r ") +e -«[2 / 8-(T'+T")] _ e -u|T'-T"|\ (3^ 



and mj? stand for the Feynman energy and Feyn- 
man mass, respectively. These definitions are under- 
stood as functional dependencies only, the optimal val- 
ues of parameters v and w are yet to be determined. 
Both Ep and mp contain terms that originate from the 
last fraction in the action l|25|l . The critical question is 
about the order of their contribution in the low temper- 
ature limit f3 — ► oo . In Eq. (|29|> the two double inte- 
grals are 0((3) and 0(1), respectively. The second one 
can therefore be omitted. The same is true about the 
two double integrals in Eq. (|3l7|) . In the second one, 
the integration region is limited to where the combi- 
nation {(3 - \t' - t"\) < (fifi)" 1 = 0(1). Since the 
pre-exponential factor is exactly the same combination 
squared, it is 0(1) throughout the essential integration 
region. As a result, the entire second double integral 
is 0(1). Therefore it can be omitted in favor of the 
first integral, which is 0(/3). Now recall that the pre- 
exponential factor (J3— \t' — t"|) 2 derives from the shifted 
boundary conditions of the phonon integration. Without 
the latter, the pre-exponential factor would have been 
|t' — t"| 2 , that is the same as in the first integral. There- 
fore, it would have been 0(/3 2 ) in the essential integra- 
tion region. That would have resulted in the second 
integral being 0(/3 2 ) and its contribution to the mass 



0(/3), which, of course, makes no sense. Thus the shifted 
boundary conditions in phonon integration are essential 
for the correct form of the pre-exponential factor and, in 
the end, for the consistency of the entire mass calculation. 
Once the second integrals in Eqs. (|29|l and l|3U|) are safely 
omitted, the first ones can be transformed according to 
the relation 

'\r'dr"g{\r'-r"\)^2^ f drg{r), (32) 
0^0 Jo 

which is valid for any function g in the limit (3 — > oo. 
The Feynman energy and mass assume their final forms: 



E h 



3 (v — w) z 



a 



(nny- 



e -hQ.T 

dr r 

„ [F(rp 



mp = mil 



^={nn)i r dr 

3V 71 " Jo 



F(r) 



1-e- 



(33) 



(34) 



(35) 
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The last step is choosing optimal values of v and w. 
The standard approach has always been to minimize Ep 
first and then substitute the obtained values into the ex- 
pression for mp. However, the variational theorem i|28|) 
is valid for any Ar and not for just the zero shift vec- 
tor. The question therefore is whether the other terms 
in the polaron action change the optimal values of the 
variational parameters. Let vq and wq minimize function 
Ep(v,w) defined by Eq. ({3"3"|> . As such, they satisfy the 
equations dE F /dv = and dE F /dw = and a neces- 
sary concavity condition. With other terms included, one 
has to minimize Ep + const//3 + (mp/2h 2 )(Ar / (3) 2 + . . . 
Note that the second and third terms are of the same 
order 0(/3 _1 ). Minimization yields v — vq + Sv, and 
w = wq + Sw, where Sv, Sw — 0(f3~ 1 ). Some parts of the 
corrections Sv and Sw depend explicitly on Ar. However, 
the minimum energy itself receives a correction that is 
only quadratic in Sv and Sw. The leading Ar correction 
to the minimal polaron action is C((Ar/ /3) 2 ), which docs 
not affect the effective mass term — (mp/2h 2 )(Ar) 2 / (3 in 
the (3 — > oo limit. 

To conclude, the polaron mass is still determined by 
the original Feynman procedure: minimize the energy 
(|33|l first, and then use the optimal parameters to calcu- 
late the mass from Eq. (|34|1 ■ 

V. SUMMARY 

In this paper, calculation of the Frohlich polaron ef- 
fective mass has been analyzed for robustness to bound- 



ary conditions in imaginary time. It has been shown 
that a consistent mass calculation must involve a corre- 
lation between the boundary conditions of the electron 
and phonon coordinates. The t = (3 end points of all 
paths have to be shifted by the same vector Ar relative 
to their r = end points. Then the effective mass is 
found as inverse diffusion coefficient of the many body 
path, cf. Eq. dJ, where Ar is used as a diffusion dis- 
tance and (3 — > oo as a diffusion time. This conclusion 
is not limited to the polaron system but is valid for any 
composite non-relativistic quantum particle. 

It has also been shown that the correlation between 
the electron and phonon boundary conditions is critical 
for the consistency of the polaron mass calculation. It re- 
sults in the correct form of a prefactor in the intermediate 
expression for mass, Eq. 1|3U|) . which allows dropping the 
second double integral altogether. In the original Feyn- 
man calculation, the second double integral was omitted 
from the beginning, i.e., during the energy calculation, 
and therefore did not to cause any problems in the mass 
calculation. Such an approach was internally inconsis- 
tent. It has been shown here that the consistency is 
restored via correct treatment of the shifted boundary 
conditions. 

Finally, it has been shown that all the above consider- 
ations do not change the numerical values of the polaron 
effective mass obtained irii^. The optimal values of the 
variational parameters are still determined by minimiza- 
tion of the ground state energy 133f) alone. The effective 
mass then follows from expression (|34|l evaluated at the 
optimal values of v and w. 
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